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We propose a set of scanning tunneling microscopy experiments in which the surface of super- 
conductor is scanned by a superconducting tip. Potential capabilities of such experimental setup are 
discussed. Most important anticipated results of such an experiment include the position-resolved 
measurement of the superconducting order parameter and the possibility to determine the nature of 
the secondary component of the order parameter at the surface. The theoretical description based 
on the tunneling Hamiltonian formalism is presented. 



I. INTRODUCTION 

Recent technological advances allowed Pan et al. to 
conduct very low temperature scanning tunneling mi- 
croscopy (STMl experiments reaching temperatures as 
low as 220 mKLl and making possible the direct imag- 
ing of the surface of high-T c superconductors (SC) in 
a superconducting state with high spatial and energy 
resolution. The Pt/Ir tips used in these experiments 
were not superconducting, but the same group has re- 
ported earlier ,0 that they were able to obtain the atomic 
resolution images of superconducting NbSe2 with the su- 
perconducting atomically sharp Nb tip in the quasipar- 
ticle tunneling regime. The idea of STM experiment in 
which the Josephson effect would be observed between 
the tip and the surface was also suggested!} 

If both the tip and the surface are superconducting 
and the temperature is low enough to have a well-defined 
phase difference between the superconductors, the tun- 
neling contact may be considered a Josephson junction 
with tunable parameters (hence the term "Josephson 
STM" or JSTM). The phase-dependent supercurrent will 
not be washed out by thermal fluctuations if the temper- 
ature is sufficiently low, 



/ c $ > 2irk B T, 



(1) 



where I c is the critical Josephson current and $o = h/2e 
is the unit magnetic flux. Assuming that the Josephson 
current can be roughly estimated as / = <rA/ep where 
A is the magnitude of the superconducting order param- 
eter (OP) and a is the normal state conductance of the 
junction, Eq. (|l|) translates into the condition imposed 
on the normal state resistance: 



R< 



A 



2k B T 



zRo 



(2) 



Here Rq = h/e 2 w 4kf2 is the resistance quantum. For 
conventional superconductors (A ~ 2 meV) and temper- 
ature T ~ 0.1 K the junction resistance should not exceed 
~ 500 KS1 for the Josephson current to be observable. For 
higher gap values, such as in high-T c materials, we esti- 
mate R < 7.5 Mil. Since the temperature enters into 
criteria (||) only via ratio A/kT, it might be easier to 
observe Josephson current using high-T c material with a 



larger value of A as a tip rather than to lower tempera- 
ture by an order of magnitude. 

Typical resistances of STM experiments are at least by 
an order of magnitude larger, than the estimated upper 
bounds for Josephson effect. Thus, Josephson supercur- 
rent is difficult to observe in a setup with an atomically 
sharp superconducting tip. Therefore we propose the 
use of a tip with a finite tunneling area. While it will 
be impossible to achieve the atomic-scale spatial resolu- 
tion, to capture the microscopic variations of supercon- 
ducting OP only the resolution comparable to the coher- 
ence length of SC is needed. This characteristic length is 
about 20 A for the high-T c materials and a few thousand 
A for conventional superconductors. Also, the finite tun- 
neling area will lead to relatively big capacitance of the 
junction, minimizing the charging energy and making the 
effects of Coulomb blockade negligible. 

Since the supercurrent depends on the magnitudes of, 
the superconducting OP on both sides of the junction,EJ 
it should be possible to measure the microscopic varia- 
tions in OP depending on the position of the tip. Such 
measurement can provide a direct insight into the micro- 
scopic nature of the superconducting state. 

Another potential application of the JSTM is the prob- 
ing of the symmetry of superconducting OP. It is known 
that a perturbation, such as an applied magnetic field 
or a surface, can induce a secondary component of the 
superconducting OP. The type of the secondary com- 
ponent can be determined from the modification of the 
Josephson current. In case when both the tip and the 
surface are unconventional superconductors, like high-T c 
cuprates, the OP is angle-dependent and the variations 
in the tunneling supercurrent may be caused by changing 
the mutual orientation, providing the direct information 
about the symmetry of OP. 

Finally, the JSTM could provide information about 
the nature of the gap in high-T c superconductors. It 
is believed, that the normal state of the underdoped 
high-r c cuprates is a pseudogap (PG) state with par- 
tially gapped Fermi surface, from which superconductiv^ 
ity emerges as the temperature is lowered through T c .El 
Currently there is no general agreement on the origin of 
PG state. Some models attribute PG to superconduct- 
ing phase fluctuations above T c p others -.-.to a compet- 
ing non-superconducting order parameter J3 Since JSTM 
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is sensitive only to the superconducting OP, the mea- 
surements of the critical Josephson current on differently 
doped samples can shed some light on the origin of PG. 
The rest of the paper is organized as follows: in Sec- 
tion |ll| we describe the tunneling Hamiltonian formalism 
used in the calculations and two practically important 
cases of JSTM, followed by proposed experimental setup 
(Section H|) and summary (Section |lV|). 



II. THEORETICAL DESCRIPTION 



We base our theoretical description of JSTM on the 
tunneling Hamiltonian formalism.Lfa In this formalism, 
time-dependent Josephson- current through the tunnel- 
ing contact is expressed aa3 



Ij(t) = 2eIm[cxp(-2ieVt/h)$ Tet (eV)}, 



(3) 



where V is the applied voltage and ^ ve ^(eV) is the 
voltage-dependent retarded correlation function, which 
in the limit of the zero temperature has the form 



*ret(^) = ^T kp T_ k ,_ 



k P 



p A k A F 

EkE p 



(4) 



surface, induces-fU-secondary component of the super- 
conducting OPEjEj Tunneling between pure s- and d- 
wave superconductors is expected to be zero by sym- 
metry arguments, at least if one pdpes not take into 
account the highep-prder processes .113 Experimental ev- 
idence exists Jl3't3ll3 that the secondary component in 
high-T c cuprates possesses s-wave symmetry In order 
to determine the dependence of the critical Josephson 
current both on the magnitude and the phase of the 
secondary OP component, we have chosen the following 
form for the superconducting OP on the surface: 



A r 



A 2 



A!, cos 2(p - 



'A 



(«) 



(5) 



Here A 2 and A), are the magnitudes of the primary 
and secondary components (respectively) and ip is the 
azimuthal angle corresponding to momentum p in the 
coordinate system associated with the surface plane. By 
varying angle a the secondary component can be assigned 
an arbitrary phase: In particular a = corresponds to 
pure d + s symmetry and a = 7r/2, to d + is. The excita- 
tion energies E k and E p are given by familiar relations 



£ k = 



g 



A? 



(6) 



and 



eV + E k + E p + iS eV - E k - E p + iS 



Here indices k and p refer to the quasiparticle momenta 
of the tip and surface, respectively, T kp is the tunneling 
matrix element, A k and A p are the superconducting or- 
der parameters, -E k and E p are the excitation energies, 
and the infinitesimal quantity iS ensures the convergence 
of the integral. 

Using Eqs. (|j) and (|J) we can calculate the Joseph- 
son current for a tunneling contact of two superconduc- 
tors. Unique feature of the Josephson effect is that below 
a certain threshold value, the dc supercurrent can flow 
without any voltage drop on the contact. This threshold 
value - critical Josephson current - is easily measured, 
because when the value of the current reaches the critical 
value, the finite voltage drop abruptly appears at the con- 
tact. If the analytical expression for the critical current 
is known, its measurements can provide the information 
about the magnitudes of OPs, its symmetry and possibly 
other quantities of interest. 



A. s-d junction 

First we consider the case when the STM tip is a con- 
ventional s-wave SC (superconducting metal, like Nb) 
with momentum-independent real order parameter A k = 
Ai and the surface under study is the a-b plane of a 
d-wave unconventional SC (e.g., YBCO or BSCCO). It 
is generally agreed, that perturbation of the unconven- 
tional d-w&ve superconductors, like magnetic field or the 
distortion of crystal lattice caused by presence of the 



£ 



(7) 



where £ k and £ p are the single-particle energies, mea- 
sured with respect to corresponding chemical potentials. 
We assume here that the tunneling matrix element is in- 
dependent of momenta, i.e. 



TkpT_k,- P = \T 



2„J<J 



(8) 



Here |T| 2 is a constant, depending only on the geometry 
of the tunneling contact (it drops exponentially with in- 
creasing distance between tip and the surface) and <p is 
the phase difference, giving rise to Josephson supercur- 
rent. 

In order to calculate the sums in Eq. (0), we convert 
them into integrals over the energies £ k , £ p and the az- 
imuthal angle (p. In the case of zero applied voltage, the 
integral with respect to energies can be calculated using 
the result due to Ambegaokar and Baratofiu to yield 



Wo) = 



2ne 2 



i<t> 



A, A, 



Ai 



tK 



Ai 



d(f 

(9) 



Here K(x) is the complete elliptic integral of the first 
kind and <7o is the conductance of the tunneling contact 
in the normal stateEl 



cr a = Aire 2 N 1 N 2 \T\ 2 



(10) 



with N\ and N 2 being the densities of states at the Fermi 
levels of tip and surface respectively. The critical Joseph- 
son current in this case is 



■fclx(7i «>«)!) 



(11) 
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where 7 = A 2 /Ai and k = A 2 /Ai are dimensionless 
parameters, 7° = Ai<j /ire and x(l> K i a ) is the value 
of integral in Eq. (||). Assuming, that Ai is fixed, we 
have studied numerically the dependence of the critical 
current on 7, k and a for realistic parameter ranges. 

The dependence of the critical Josephson current for 
a = (corresponds to pure d + s symmetry) on 7 and k 
is shown on Fig. [II 
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FIG. 1. The dependence of the critical Josephson current 
on 7 and k for an s-d junction in the limit of zero applied 
voltage and a = 0. Different curves correspond to different 
values of k, starting from k = 0.01 (lowest curve) with the 
increment of 0.01. The uppermost curve thus corresponds to 
K = 0.1. 



Another important case is a — 7r/2, which corresponds 
to d + is symmetry of superconducting OP. The results 
for this case are presented in Fig. 0. 
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FIG. 2. The dependence of the critical Josephson current 
on 7 and k for an s-d junction in the limit of zero applied 
voltage and a — n/2. Different curves correspond to different 
values of k, starting from k — 0.01 (lowest curve) with the 
increment of 0.01. The uppermost curve thus corresponds to 
« = 0.1. 

Useful information can also be extracted from the ob- 
servation of quasiparticle current between superconduc- 
toxs_at a finite temperature. The current can be written 
asfflJ 



2?re 



dip 



+ OC 



\E\ 



\E + eV\ 



^E* - A\ ^{E + eVf - |A 2 



[f(E)-f(E + eV)]dE, 



(12) 



where f(E) is a Fermi-Dirac distribution function and 
the regions \E\ < Ai and \E + eV\ < |A 2 | are excluded 
from the integration. We have numerically calculated 
current- voltage characteristic for a realistic set of param- 



15 



W/aW 



0.025. 



eters: Ai/fcTj= 2, A 2 d) /Ai = 10 and A 2 S 7A 2 
It is knownJ13 that when the tunneling contact of two 
s-wave superconductors with different constant real or- 
der parameters Ai and A2 is considered, the quasipar- 
ticle current-voltage characteristic possesses two signifi- 
cant features: a peak at the voltage |Ai — A 2 | and an 
abrupt increase in current at voltage Ai + A 2 , corre- 
sponding to transition to Ohmic regime. We find that the 
position of the peak and the Ohmic transition depends on 
the symmetry of the superconducting OP. In particular 
(as seen in Fig. 0) the current is sharply peaked at the 
voltage A 2 — Ai for a = n/2 (d+is symmetry), while 

for a — the peak is at the voltage A 2 + A 2 — Ai. The 
transition to Ohmic regime is also shifted to higher volt- 

( s) 

ages by the value of A 2 . Thus, such a measurement can 
complement the measurements of Josephson critical cur- 
rent to provide additional information about the nature 
and magnitude of superconducting OP on the surface. 
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FIG. 3. Typical quasiparticle current-voltage characteristic 
(values of parameters are given in the text) of a tunneling con- 
tact between d-wave superconductor with the secondary OP 
component and s-wave superconductor. Different values of a 
correspond to d+s (a — 0) and d+is (a = tt/2) symmetries. 
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B. d-d junction 

The d-wave nature of a pair requires finite size tunnel- 
ing area. If the tunneling occurs from only one site on 
the tip it would imply effective momentum averaging of 
pair wavefunction and would yield zero Josephson cur- 
rent for d-wave pairs. Thus it is desirable to have a finite 
region of the tip where tunneling can occur. In this case, 
one looses spatial resolution, compared to conventional 
tips. However the gain is the phase coherent current be- 
tween the tip and the surface. While the idea of using 
the high-T c cuprates as a material of STM tip has not 
yet been reported in the literature, the fabrication of rel- 
atively small (of the order 100 A) flat high-T c probes, 
suitable for use as tips, is feasible. 

Theory based on tunneling Hamiltonian predicts, 
that if the tunneling matrix element is momentum- 
independent, c-axis Josephson current at zero temper- 
ature is identically zero for superconductors with pure 
rf-wave symmetry of OP. It is possible to argue that, as 
discussed in previous section, the presence of the surface 
will induce the secondary OP component, which can lead 
to a nonzero current. However, since the magnitude of 
secondary component is usually much smaller than the 
magnitude of primary one, this induction will produce 
a second order effect and we assume that it can be ne- 
glected. 



hypothesis leads to non-zero current even in the case of 
pure d-wave symmetry. 

We consider an a-b plane tunneling contact of two d- 
wave superconductors. We allow the magnitudes of the 
superconducting OP to be different. Because of the angu- 
lar dependence of the superconducting OP, the Josephson 
current will depend on the mutual orientation of two su- 
perconductors and we introduce angle 9 (angle between 
a-axes of superconductors) to describe this orientation. 
The problem is essentially two-dimensional and the or- 
der parameters of the tip Ai and surface A2 may be 
written as 



and 



Ai = A^ 0) cos2<^ 



A 2 = A^ 0) cos2(^ + 6»). 



Tunneling matrix element has the form 
T kp T_ k) _ p = |T|V^(k-p), 



(13) 



(14) 



(15) 



where k and p are two-dimensional vectors. Again, we 
change the sums in Eq. (|j) into integrals, which can be 
calculated numerically. In the case of zero applied volt- 
age and zero rotation angle 9 it is possible to calculate the 
integral analytically and find that the critical Josephson 
current in this limit is 




Ir. = I, 



7hi7 2 



1' 



(16) 



where 7 = A^ 0) /A^ 0) , 



a(°) 
assuming l\\ 



M 



< A ( "> and J° is 



constant current independent of 7. This dependence is 
presented in Fig. |[ 7 = corresponds to absence of the 
superconductivity in the tip, which leads to dissapear- 
ance of the Josephson current. 

We have also studied the dependence of the critical 
current on the rotation angle 9 at zero applied voltage 
(Fig. g). As expected, it has a maximum at 9 = 0, 
monotonically decreases with increasing 9 and becomes 
zero at 9 = ir/4. For the angles bigger than 7r/4, the 
critical current can be obtained as 



FIG. 4. The dependence of the critical Josephson current 
on the ratio of OP magnitudes for a d-d junction in the limit 
of zero applied voltage and zero rotation. 



Here we discuss a different scenario, supported by 
recent—theoretical and experimental work of Latyshev 
et alxJl In their work, the intrinsic c-axis quasiparti- 
cle transport between layers of BSCCO was studied. 
It was found that interlayer transport is predominantly 
momentum-conserving, i.e. the in-plane components of 
quasiparticle momentum are conserved in the tunneling 
process. It seems natural to make an assumption, that 
in-plane momentum is also conserved in case of Joseph- 
son supercurrent. As we show below, adoption of this 



'<G 



= -*■' — ' 



(17) 



i.e. for 9 > 7r/4 the current changes direction, but 
the curve remains symmetric with respect to the point 
9 = 7r/4. It is notable, that since the curve is symmet- 
ric with respect to 9 = 0, there is a kink in the angle 
dependence at 9 = 0. 

It can be seen from Fig. g, that the decrease is almost 
linear with 9 and linear approximation becomes better 
as 7 is reduced. For practical purposes an approximate 
formula 



Ir = l 



o7 m 7 



4 
,1 < 

IV 7T 



(18) 



may be used. 
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FIG. 5. The dependence of the critical Josephson current 
on the rotation angle 9 for a d-d junction in the limit of zero 
applied voltage at different values of j (solid lines). Dashed 
lines show the linear approximation of Eq. (UM . 



It is worth mentioning that when the angle 9 is non- 
zero, the denominator in Eq. (jj) has a positive lower 
bound. This means that for voltages much smaller than 
this lower bound, the integrand may be expanded in pow- 
ers of voltage. It turns out that all first order terms in 
voltage exactly cancel out and dependence on voltage 
enters only as higher order corrections (proportional to 
the square of the voltage), so the zero- voltage behavior 
is not altered significantly. The range of applicability of 
this approximation depends on actual values of 7 and 9. 



III. EXPERIMENTAL SETUP 

The sketches of experimental setup for s-d and d-d tun- 
neling are shown in Fig. rc| and Fig. ffl respectively. In 
both cases, superconducting tip and surface form an elec- 
tric circuit together with current (voltage) source and 
measuring equipment. 




the normal state conductance of the tunneling junction 
at the specific position on the surface is measured. Next, 
the measurement of the critical Josephson current (by 
starting with zero current and increasing it until a volt- 
age drop appears on the contact) is performed. Using 
this information, a conclusion can be made about the 
magnitudes of the primary and/or secondary OP com- 
ponents on the surface. By repeating the measurement, 
a surface map of variations in order parameter may be 
obtained. Both the surface, formed by a-b plane and the 
one, parallel to c-axis can be studied. 




FIG. 6. Proposed experimental setup for s-d tunneling. 



FIG. 7. Proposed experimental setup for d-d tunneling. 



In d-d case, the apparatus must be designed in such 
a way that the tip can be rotated around vertical axis. 
Measuring the Josephson current as a function of rotation 
angle, the information about the symmetry of OP can be 
extracted. If Josephson current has nodes, separated by 
A9 = 7r/2, a conclusion can be made that the non-zero 
Josephson current is due to coherent tunneling and su- 
perconductors possess pure <i-wave symmetry. If, on the 
other hand, current never becomes zero or becomes zero 
with different periodicity, then the secondary component 
of the OP must be a primary reason for the existence of 
the current. 

Finally, the experiments aimed at establishing the na- 
ture of the pseudogap should be conducted with differ- 
ently doped tip and surface. For example, measuring 
critical Josephson current in a setup with underdoped 
tip and overdoped surface and comparing it with mea- 
surements on identically underdoped (or overdoped) con- 
tacts will provide information about the magnitude of the 
pseudogap and whether or not it is of superconducting 
origin. Detailed discussion of such an experiment will be 
given elsewhere. 



IV. SUMMARY 



In s-d case, measurements can be performed in two 
stages. First, by driving a large current through the tun- 
neling contact, the superconductivity is suppressed and 



In summary, we have a proposed a series of STM exper- 
iments of a superconducting surface with a superconduct- 
ing tip. Treating the tunneling contact as a Josephson 
junction with tunable parameters, we have developed a 
theoretical description of these experiments, calculating 
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the dependence of critical Josephson current on various 
relevant parameters. Since Josephson current is sensitive 
to the magnitudes of superconducting order parameters 
in the tip and in the surface, this new technique allows to 
measure the microscopic spatial variations of the OP on 
the surface and determine the magnitude and phase of 
the secondary OP component, induced by the surface. It 
is also possible to probe the symmetry of superconduct- 
ing OP directly, exploiting the dependence of the Joseph- 
son critical current on the mutual orientation of rf-wave 
superconductors. Finally, since the Josephson current is 
sensitive only to the superconducting gap, the new tech- 
nique can provide information about the nature of the 
pseudogap, in particular, confirm or falsify its supercon- 
ducting origin. 
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